Relativistic gravitational deflection of light and its impact on the modeling accuracy 

for the Space Interferometry Mission 
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We study the impact of relativistic gravitational deflection of light on the accuracy of future 
Space Interferometry Mission (SIM). We estimate the deflection angles caused by the monopole, 
quadrupole and octupole components of gravitational fields for a number of celestial bodies in 
the solar system. We observe that, in many cases, the magnitude of the corresponding effects is 
significantly larger than the 1 ^as accuracy expected from SIM. This fact argues for the development 
of a relativistic observational model for the mission that would account for the influence of both 
static and time-varying effects of gravity on light propagation. Results presented here are different 
from the ones obtained elsewhere by the fact that we specifically account for the differential nature 
of the future SIM astrometric measurements. We also obtain an estimate for the accuracy of possible 
determination of the Eddington's parameter 7 via SIM global astrometric campaign; we conclude 
Oh' that accuracy of ~ 7 x 10^® is achievable via measurements of deflection of light by solar gravity. 
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I. INTRODUCTION 

'o' 

r-jH' 

y ' The last quarter of the 20th century has changed the status of Einstein's general theory of relativity from a purely 
^ . theoretical discipline to a practically important science. Today general relativity is the standard theory of gravit 
especially where the needs of astronomy, astrophysics, cosmology and fundamental physics are concerned [if, 0, 
As such, this theory is used for many practical purposes involving spacecraft navigation, geodesy and time transfer. 
T-H ■ Present accuracy of astronomical observations already requires relativistic description of light propagation as well as 
K*" [ the relativistically correct treatment of the dynamics of the extended celestial bodies As a result, some of the 
leading static-field post-Newtonian perturbations in the dynamics of the planets, the Moon and artificial satellites 
have been included in the equations of motion, and in time and position transformation. It is also well understood 
that effects due to non-stationary behavior of the solar system gravitational field as well as its deviation from spherical 
symmetry should be also considered [51] and implemented in the appropriate models. 

Space-based astrometry has brought about a renaissance in the entire field of astrometry that is perhaps the most 
fundamental, and oldest of all areas in astronomy The ESA Hipparcos mission, which operated from 1989-1993, 
yielded an astrometric catalog of 118,000 stars down to 12.5 magnitude, with positional accuracy of 1 mas for stars 
brighter than V = 11. The European Space Agency (ESA) is now developing the Gaia mission as a next generation 
astrometric survey mission 0, Q, which is expected to produce a catalog of ^ 10^ stars, with accuracy ~ 20-25 
microarcsec (fias) for stars brighter than V = 15. Precision astrometry remains a cornerstone of the field and is 
' poised to make a major impact on many fields of modern astronomy, astrophysics, and cosmology 
5^ , NASA's SIM PlanetQuest mission, hereinafter SIM, as another example of a space-based facility instrument for 
astrometry. The acronym SIM stands for Space Interferometry Mission. SIM will be the first space-based Michelson 
interferometer for astrometry. The instrument will operate in the optical waveband using a 9-m baseline between the 
apertures. With a global astrometry accuracy of 3 /las for stars brighter than V = 20, it will measure parallaxes and 
proper motions of stars throughout the Galaxy with unprecedented accuracy. Operating in a narrow-angle mode, it 
will achieve a positional accuracy of 0.6 /las for a single measurement, equivalent to a differential positional accuracy 
at the end of the nominal 5-year mission of = 0.1 /las. This performance is about 1000 times better than existing 
capabilities on the ground or in space, and about 100 times better than the upcoming Gaia mission, for differential 
measurements. Such high accuracy will allow SIM to detect and measure masses of terrestrial planets around stars 
in our Galactic neighborhood (see [6l| for review). 

SIM is a targeted mission which measures the astrometric positions of stars, referencing the measurements to a grid 
of 1302 stars covering the entire sky. Its scheduling is highly flexible, in both the order of observations, their cadence, 
and the accuracy of each individual measurement. This contrasts with the Hipparcos and Gaia missions, which scan 



* Email: turyshev@jpl.nasa.gov, turyshev@sai.insu.ru 



2 



Apparent position of 

tlie reference source 
.r 



Reference 



the entire sky according to a pre-determined scanning pattern. Many astrometry experiments can make effective use, 
or in some cases require, this pointing capabiUty for instance, searches for terrestrial planets (especially in multiple 
planet systems), stellar microlensing events, orbits of eccentric binary systems, and variable targets such as X-ray 
binaries and active galactic nuclei. Currently, the ICRF, defined by the locations of 212 extragalactic radio sources 
[1) [ifl] with most having errors less than 1 mas, is the standard frame for astrometry. SIM is expected to yield an 
optical reference frame at a level of about 3 /ias; it will be tied' to the ICRF by observing a number of radio-loud 
quasars in common. 

In this paper, we discuss the rapidly 
forming field of applied general relativity 
to demonstrate its influence on the high- 
precision astrometry. Recent advances in 
the accuracy of astrometric observations have 
demonstrated importance of taking into ac- 
count the relativistic effects introduced by the 
solar system's gravitational environment. It 
is known that the reduction of the Hipparcos 
data has necessitated the inclusion of stellar 
aberration up to the terms of the second or- 
der in v/c, and the general relativistic treat- 
ment of light ben ding due to the gravitational 
field of the Sun [uf and Earth Even 
higher modeling accuracy is anticipated for 
Gaia [13]. 

Prediction of the gravitational deflection of 
light was one of the first successes of gen- 
eral relativity. Since the first confirmation by 
the Eddington's expedition in 1919, the effect 
of gravitational deflection has been studied 

quite extensively and currently analysis of al- F^G. 1: Geometry of gravitational deflection of starlight by the Sun. 
most every precise astronomical measurement must take this effect into account [l|, [3j [Hi [3 ■ According to general 
relativity, the light rays propagating near a gravitating body are achromatically deflected by the body's relativistic 
gravity field. The whole trajectory of the light ray is bent towards the body by an angle depending on the strength 
of the body's gravity. The solar gravity field produces the largest effect on the light traversing the solar system. To 
first order in the gravitational constant, G, the solar deflection angle 0®. depends only on the solar mass Mq and the 
impact parameter d relative to the Sun flT'l : 
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= 1.751 



d 



1 + cos X 



(1) 



where Rq solar radus. The absolute magnitude for the light deflection angle is maximal for the rays grazing the sun, 
e.g. 0®. = 1.751 seconds of arc. Most of the measurements of the gravitational defiection to date involved the solar 
gravity field, planets in the solar system or gravitational lenses. For the future astrometric observations with SIM, in 
addition to the Sun, effect of planetary gravitatiobal defiections of light must also be considered. The most precise 
measurement of the light defiection with the planet Jupiter was done in |l^ . Relativistic deflection of light has been 
observed, with various degrees of precision, on distance scales of 10^ to 10^^ m, and on mass scales from 10~^ to lO'^'^ 
solar masses, the upper ranges determined from the gravitational lensing of quasars [l9l.[20j. 

In the case of SIM, the star is assumed to be at a very large distance compared to the Sun, and x is the angular 
separation between the deflector and the star. With the space observations carried out by SIM, x is not necessarily 
a small angle. The relevant geometry and notations are shown in Fig. [T] In this figure we emphasized the fact that 
the difference of the apparent position of the source from it's true position depends on the impact parameter of the 
incoming light with respect to the deflector. For the astrometric accuracy of a few fias and, in the case when the Sun 
is the deflector, positions of all observed sources experiencing such a displacement. This is why, in order to correctly 
account for the effect of gravitational defiection, it is important to process together the data taken with the different 
separation angles from the defiector. In the wide-angle astrometry mode SIM, will be observing the sky in a 15° 
patches of sky (called field of regard or FoR) making a set of differential observations within the FoR. Therefore, this 
differential nature of the measurements would result in minimizing the contribution of the gravitational defiection on 
the single measurement. To refiect this fact, we will present results for two types of astrometric measurements, namely 
for the absolute (single ray deflection) and differential (two sources separated by the 15° field of regard) observations. 

A major objective of this paper is to show that, before microarcsecond-level astrometry will become a powerful 



3 



tool for 21-st century astronomy, there is a need for an adequate modeling necessary to match this new frontier of 
astrometric accuracy. The prospect of new high precision astrometric measurements from space with SIM requires 
inclusion of relativistic effects at the {v/cY level [a,l2l|. At the level of accuracy expected from SIM, even more subtle 
gravitational effects on astrometry from within the solar system will start to become apparent, such as the monopole 
and the quadrupolc components of the gravitational fields of the planets (2^ and the gravito-magnetic effects caused 
by their motions and rotations. Thus, the identification of all possible sources of "astrophysical" noise that may 
contribute to the future SIM astrometric campaign, is well justified. 

This work is organized as follows: Section |TT] discusses the influence of the relativistic deflection of light by the 
monopole components of the gravitational fields of the solar system's bodies. We present the model and our estimates 
for the most important effects that will be influencing astrometric observations of a few /^as accuracy, that will be 
made from within the solar system. Section Hill will specifically address three most intense gravitational environments 
in the solar system, namely the vicinities of the Sun, Jupiter and Earth. In Section [IVI we will discuss the effects of 
the gravitational deflection of light by the higher gravitational multipoles (both mass and current ones) of some of 
the bodies in the solar system. We derive constraints on the navigation of the spacecraft and the accuracy of the 
solar system ephemerides. In Section |V] we investigate the possibility of improving the accuracy of the Eddington's 
parameter 7 via astrometric tests of general relativity in the solar system. We also discuss the opportunity to measure 
the solar acceleration towards the Galactic Center with SIM. We will conclude the paper with the discussion of the 
results obtained and our recommendations for future studies. 



II. GRAVITY CONTRIBUTIONS TO THE LOCAL ASTROMETRIC ENVIRONMENT 

In this Section we develop a model for light propagation that will be used to estimate various relativistic effects 
due to gravitational deflection of light by the solar system's bodies. 



A. Relativistic deflection of light by the gravity monopole 



The first step into a relativistic modeling of a light path consists of determining the direction of the incoming 
photon as measured by an observer located in the solar system as a function of the barycentric coordinate position of 
the light source. Apart from second and third orders of velocity aberration the only other sizable effect is due to the 
bending of light rays in the gravitational field of solar system bodies [2l|, [2^ . Effects of the gravitational monopole 
deflection of light are the largest among those in the solar system. 

Generalizing on a phenomenological parameterization of the gravitational metric tensor field, which Eddington 
originally developed for a special case, a method called the parameterized post-Newtonian (PPN) formalism has been 
developed (see [4I for discussion). This method represents the gravity tensor's potentials for slowly moving bodies 
and weak inter-body gravity, and is valid for a broad class of metric theories, including general relativity as a unique 
case. The several parameters in the PPN metric expansion vary from theory to theory, and they are individually 
associated with various symmetries and invariance properties of the underlying theory (see [l3l for details). 

If (for the sake of simplicity) one assumes that Lorentz invariance, local position invariance and total momentum 
conservation hold, the metric tensor for a single, slowly-rotating gravitational source is given by: 

„GA/ / ^ i?2 3cos2 6l- 1\ 4, 

go. = 2(7 + l)^i|^^+0(c-^), (2) 
5., = -S.,[l + 2j—[l-.h — )+0(c 5), 

where M and S being the mass and angular momentum of the body, J2 and R are the body's quadrupolc moment and 
its radius, and G is the universal gravitational constant, r is the distance from the center of the body to a particular 
point. The term in 500 is the Newtonian limit; the 1/c^ term in goi and the term in gij, are post-Newtonian 
corrections. All of these terms are required to describe light propagation phenomena to the first post-Newtonian 
order. 

The Eddington parameter 7 in the Eqs. ([2]) represents the measure of the curvature of the space created by a unit 
rest mass [ITj. Note that general relativity, when analyzed in standard gauge of the PPN formalism (see 
for details), gives: 7=1. The Brans-Dicke theory is the most famous among the alternative theories of gravity. 
It contains, besides the metric tensor, a scalar field 4> and an arbitrary coupling constant cu, related to this PPN 
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parameter as 7 = The stringent observational bound resulting from the 2003 experiment with the Cassini 

spacecraft require that 40000 [2^. There exist additional alternative theories that provide guidance for 

gravitational experiments |24l |. 

Metric tensor Eq. ^ can be used to derive expressions need to describe propagation of electro-magnetic signals 
between any of the two points in space. Following the standard procedure of integrating light geodesies (see [3, HB] 
for details), the corresponding light-time equation for a single deflecting body can be derived in the following form 



^ + (1 

c ' c 



GM 
7)^^ In 



ri+r2 + ri2 



ri + r2 - ri2 



(3) 



where ti refers to the signal transmission time, and t2 refers to the reception time. ri^2 a-re the barycentric positions 
of the transmitter and receiver, and ri2 is their spatial separation (see [2^ for details). Also, for the moment, we 
neglected the presence of J2 and S in the Eqs. but will investigate contributions of these and other multipoles in 
SecHVl 

For a realistic observing scenario with the SIM, the sources of light are located far out side the solar system, 
r2 ri = rs and ri2 can be approximated as ri2 = jri — r2| ~ ri — {ns ■ 7^2), where we introduced a notation 
fis = fi /ri . This approximation allows one to represent the expression in the square brackets in Eq. ([3]) as follows: 



ri+r2 + ri2 ^ 2ri + r2 - (ng ■ f2) 
r2 + {ns ■ ^2) 



(4) 



External Path Delay 
d = Bcos0 + C 



telescope 2 



ri +r2 ~ ri2 

The SIM instrument does not directly measure the 
angular separation between stars, but the projection 
of each star direction vector onto the interferometer 
baseline by measuring the pathlength delay of starlight 
as it passes through the two arms of the interferome- 
ter. The SIM instrument will precisely measure optical 
path difference (OPD) between the wavefronts of light 
received by the two telescopes forming the interfero- 
metric baseline (see Fig.[2]). (In fact, for SIM, with its 
10 m baseline, such a measurement will be done with 
a precision at the level of 1 picometer.) This differ- 
ence will result in the different internal OPDs needed 
to apply in order to coherently add the signals. The 
delay measurement is made by a combination of inter- 
nal metrology measurements to determine the distance 
the starlight travels through each arm, external metrol- 
ogy measurements that determine the length and lo- 
cal orientation of the baseline, and a measurement of 
the central white light fringe to determine the point 
of equal optical pathlength [gj. Therefore, the OPD 
is the main observable that the interferometric instru- 
ment will measure; relativistic modeling this delay will be among the main objectives of the upcoming SIM modeling 
effort. 

An interferometer measures optical path difference (OPD) between the wavefronts of light received by the two 
telescopes forming the interferometric baseline b at points r2 and = 7^2 + To account for that fact we need to 
determine the temporal difference between the signals received at these telescopes which is £ = c(<2 — ^i) — c{t2 — ti) = 
c{t'2 — ^2)- The first term in Eq. ([3]) is the geometric delay. Using approximations (6 ^ r2 ^ ri), it is easy to see that 
this term leads to the approximate expression for geometric delay, given as 




Internal 
Path Delay 

Delay lit 

FIG. 2: Basic geometry of light propagation in the stellar inter- 
ferometry. 
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^2-^12 = -{h-ns)/{l~ns-r2/c). 
In this paper, we concern with only the largest contributions from the gravitational defection of light, thus most of 
the velocity-dependent terms will discarded (see [H, H^l for details). 

The second term Eq. ([3]) is the relativistic delay £gr, the focus of this work. Using Eq. ([4]) one can present expression 
in the square brackets of Eq. ^ as below 



2ri + r2' - {ns ■ f^) r2 + {ns ■ ^2) 



2ri + r2 - {ns ■ r2) ^2' -(- {ns ■ rj) 



1 

1 



h{ns + ^2) 
7'2 1 + {ns ■ "2)' 



(5) 



In the first order in gravitational constant, one can add individual interferometric delays due to the gravity of the 



bodies along the light path. As a result, the general relativistic contribution to the OPD 



-gr 



CTgr takes the following 
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approximate form: 



-(7 + 1) 



GMf 



B 



TB 



h{ns + ub) 



l + {ns ■ Ub) 



(6) 



& = 6 (cose, sine), = (cosas, sinas), 715 = (cosi 



where tb is the distance from SIM to a deflecting body i?, hb = ^bI^b is the unit vector in this direction. This OPD 
is the leading general relativistic observable that the interferometric instrument will measure; a complete relativistic 
modeling this delay should be the main objective of the upcoming SIM modeling effort. 

In general, a three-dimensional approach must be used in order to work out a practical model of the interferometric 
time delay. Howevere, for the purpose of this paper, it is sufficient to confine our analysis to a plane and parameterize 
the quantities involved as follows (see Fig. [3]): 

sin0), (7) 

where e is the angle of the baseline's orientation with respect to the instantaneous body-centric coordinate frame, as 
is the right assention angle of the interferometer as seen from the this frame and B is the direction to the observed 
source correspondingly. The geometry of the problem and notations are presented in the Figure [3) 

It is convenient to express the gravita- 
tional contribution to the total OPD Eq.® 
in terms of the deflector and the source sep- 
aration angle xb as observed by the inter- 
ferometer. In our approximation, the fol- 
lowing relations d = rs sinxs = sinxs 
and xs + Xb + 9 — a — tt are valid; this 
allows one to eliminate angle xs by ex- 
pressing the source's position angle 9 via 
the separation angle xb as below: 

• r i^B 

tt+ub-Xb- 



Interferometer 



rs 



sm XB 



(8) 




Deflecting Body 



As the sources will be located at a very 
large distance, rs, compares to the distance 
between the interferometer and the deflec- 
tor {tb <C rs), we can neglect the presence 
of the last term in the equation Eq.®, so 

that ~ TT -I- ttB — Xb- After substituting expressions ([5]) and ([7]) into Eq. we rewrite the contribution of the 
gravitational deflection to the total OPD, Eq.®, in the following form: 



FIG. 3: Geometry and notations for tlie gravitational deflection of light. 



= -(7 + l)E 



G MBb 



rB 



cos(e — as) + sin(e — cub 



1 +C0SXB 



sm XB 



(9) 



This expression describes gravitational delay as measured by an interferometer; we will use it for estimation purposes. 



1. Absolute Astrometric Measurements 



Eq. ([9]) is appropriate for estimation the magnitudes of the gravitational bending effects measured by an interfer- 
ometer. It depends on the angle between the baseline and deflector-instrument vectors, e — a^. As the main objective 
of this paper is to estimate the magnitudes of the effects involved, we choose e — as = ^ that maximizes contribution 
of each individual deflector for a particular orbital position of the spacecraft and the baseline orientation. As a result, 
in the SIM proper reference frame Eq. ([9]) may be re-written as £gr = — X^b ^gri with the individual contributions of 
the deflecting bodies to gravitational delay and deflection angle 0^ ~ in the following form 

nB ( , ..G MBb 1 + COSXB , r,B / , . G Afs 1 + COS xs ,,„^ 

4r = -(7 + l)^ : and 9 =-["f+l) — ■. . (10) 

^ &■ rB smxs rs sm^B 

The two expressions £^ and will be used interchangeably throughout the paper. 

For complete analysis of the gravitational deflection of light we will have to account for the time dependency in all 
the quantities involved. Thus, one will have to use the knowledge of the position of the spacecraft in the solar system's 
barycentric reference frame, the instrument's orientation in the proper coordinate frame psj . the time that was spent 
in a particular orientation, the history of all the maneuvers and re-pointings of the instrument, etc. These issues are 
closely related to the principles of the operational mode of the instrument that is currently still being developed. 
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TABLE I: Relativistic monopole deflection of by the solar system bodies at the SIM's location. 



Solar 


Angular size 


Deflection of grazing rays 


system's 


at SIM pos., 


absolute 


Qin. iio 1 


rllflr 11"! 
Qlii. 11 1 


object 


/v£j, arcsec 


nB 

Pgr , /ias 


r/i-B 

oPgr, /.las 


oWg]-, /las 


Sun 


0°. 26656 


1". 75064 


1". 72025 


1".38221 


Sun at 45° 


45° 


9 831.39 


2 777.97 


237.66 


Moon 


47.92690 


25.91 


25.87 


25.56 


Mercury 


5.48682 


82.93 


82.92 


82.81 


Venus 


30.15040 


492.97 


492.69 


488.88 


Earth 


175.88401 


573.75 


571.90 


547.03 


Mars 


8.93571 


115.85 


115.83 


115.57 


Jupiter 


23.23850 


16419.61 


16412.60 


16 314.30 


Jupiter at 30" 


30.0 


12 719.12 


12 712.03 


12 614.21 


Saturn 


9.64159 


5 805.31 


5 804.27 


5 789.79 


Uranus 


1.86211 


2171.38 


2171.30 


2170.26 


Neptune 


1.18527 


2 500.35 


2 500.29 


2 499.52 


Pluto 


0.11478 


2.82 


2.82 


2.82 



2. Differential Astrometric Measurements 



SIM will perform its astrometric campaign working in differential mode either within FoR=15° for the wide angle 
astrometry or within FoR=l° for the narrow angle astrometry. To evaluate the impact of gravitational delay of light 
on these measurements, we need to derive the appropriate expressions reflecting the differential nature of astrometric 
measurements with SIM. 

Within the accepted approximation, the necessary expression for the differential OPD may be obtained by sub- 
tracting OPDs for the different sources one from one another. Using Eq. this results in the following expression: 



= i?,r ~ = -(7 + 1) E - ^^i^^i , (11) 



where 7151 and 7152 are the barycentric positions of the primary and the secondary objects. By using parameterization 
for the quantities involved similar to that above {b <^ rs TS1.S2), this expression may be presented in terms of the 
deflector-source separation angles, xiB, X^b, as follows: 

(^^gr = -(7 + 1) >^ — sm(e - ub) — . (12) 

^ rB sm 2XS1 sm ^xb2 

Similar to the discussion of absolute defection angles, we choose e — = ^ that maximizes contribution of each 
individual deflector for a particular orbital position of the spacecraft and the baseline orientation. Therefore, in the 
SIM proper reference frame Eq. (|12p may be re- written as Slgr — — X^b ^^gr^ with the individual contributions for 
gravitational delay <5£^ and corresponding deflection angle 69^^ ~ Si^^/b in the following form 

xeB I , G Meb sin ^ (x2B - Xib) „ b , , G Mb sm^\{X2B ~ Xib) , . 

'^^gr = (7 + l)— ^ — and (56*^ (7 + 1)— ^ — . (13) 

^ TB sm 2X1S • sm 2X2B tb sm ^xiB ■ sm ^X2B 

The two expressions Sig^. and Sdg^ will be used interchangeably throughout the paper. 



B. Deflection of Grazing Rays by the Bodies of the Solar System 



We are now ready to evaluate the influence of the solar system's gravity field on the future high-accuracy astrometric 
observations. In particular, we estimate the magnitudes of the angles of gravitational deflection for those light rays 
that are grazing the surfaces of celestial bodies. 

Table U shows magnitude of the angles characterizing relativistic monopole deflection of grazing (e^. xiB = T^b) 
light rays by the bodies of the solar system at the SIM's location (i.e., the solar Earth-trailing orbit Q)- Results for 
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TABLE II: Relativistic monopole deflection of light: the angles and the critical distances for A6o = 1 fias astrometric accuracy. 
Solar defection of light must be always taken into account, as at the SIM's position at 1 AU from the Sun, the solar gravitational 
deflection effects are always larger than 1 ^as. The critical distances for the Earth are given for two distances, namely for 0.05 
AU (27° .49) and 0.01 AU (78°. 54). 



Object 


6", aas 


Critical distances for accuracy of 1 /las 


"crit, km 


"crit, deg 


"crit , K-B 


oun 


1 . ( 0UD4 


always 


always 


always 


Moon 


25.91 


4.501 X 10^ 


0°.34 - 1°.72 


25.9 ■ Tim 


Mercury 


82.93 


2.023 X 10^ 


0°.06 - 0°.13 


82.9 • Umb 


Venus 


492.97 


2.982 X 10*^ 


0°.66 - 4°. 13 


492.9 ■ TZv 


Earth 


573.75 


3.453 X 10^ 


27° .49 - 78°. 54 


541.4 ■ 7^® 


Mars 


115.85 


3.931 X 10^ 


0°.06 - 0°.29 


115.9 ■ TZms 


Jupiter 


16 419.61 


6.270 X 10* 


64°.06 - 88°.51 


8 849 • TZj 


Saturn 


5 805.31 


3.420 X 10* 


12°. 56 - 15° .45 


5 700 ■ TZs 


Uranus 


2171.38 


5.319 X 10^ 


1°.01 - 1°.12 


2 171 • TZu 


Neptune 


2 500.35 


6.276 X 10^ 


0°.77 - 0°.82 


2 500- 71 jv 


Pluto 


2.82 


9.025 X 10^ 


0".31 - 0".32 


2.8•7^p 



absolute deflection angles agree vifith values obtained by other authors (for instance, [23j). The resuhs presented in 
the terms of the following quantities: 

i) . for absolute astrometry results are given in terms of the absolute measurements £^ and 0^ from Eq. (fTO|) : 

ii) . for differential astrometry resuhs are given in terms of the absolute measurements S£^^ and 69^^. from Eq. (jl3p . 

For the differential observations the two stars are assumed to be separated by the size of the instrument's field of 
regard. For the grazing rays, position of the primary star is assumed to be on the limb of the deflector. Moreover, 
results are given for the smallest distances from SIM to the bodies (e.g. when the gravitational deflection effect is 
largest). For the Earth-Moon system we took the SIM's position at the end of the first half of the first year mission 
at the distance of 0.05 AU from the Earth. Presented in the right column of Table U are the magnitudes of the body's 
individual contributions to the gravitational delay of light at the SIM's location. 

Note, that the angular separation of the secondary star will always be taken larger than that for the primary. It 
is convenient to study the case of the most distant available separations of the sources. In the case of SIM, this is 
the size of the field of regard (FoR). Thus for the wide-angle astrometry the size of FoR will be 15° = rad, thus 
X2B = Xi-B + For the narrow-angle observations this size is FoR = 1° = rad, thus for this type of astrometric 
observations we will use X2B = Xis + xfo! finally, & = 10 m is the baseline length used in the estimates. 



1. Critical Impact Parameter for High Accuracy Astrometry 

The estimates, presented in the Table U have demonstrated that it is very important to correctly model and account 
for gravitational influence of the bodies of the solar system. Depending on the impact parameter ds (or planet- 
source separation angle, xb), one will have to account for the post-Newtonian deflection of light by a particular 
planet. Most important is that one will have to permanently monitor the presence of some of the bodies of the solar 
system during all astrometric observations, independently on the position of the spacecraft in it's solar orbit and the 
observing direction. The bodies that introduce the biggest astrometric inhomogeneity are the Sun, Jupiter and the 
Earth (especially at the beginning of the mission, when the spacecraft is in the Earth' immediate proximity) . 

Let us introduce a measure of such a gravitational inhomogeneity due to a particular body in the solar system. To 
do this, suppose that future astrometric experiments with SIM will be capable to measure astrometric parameters 
with accuracy of A^g = M^-s, where Ak is some number characterizing the accuracy of the instrument (e.g. for 
a single measurement accuracy Afc — 3 for stars brighter than V = 20 and for the mission accuracy Ak = 0.1, see 
0). Then, there will be a critical distance from the body, beginning from which, it is important to account for the 
presence of the body's gravity in the vicinity of the observed part of the sky. We call this distance - critical impact 
parameter, df^^^, the closest distance between the body and the light ray that is gravitationally deflected to the angle 

eUdgit) = A^o = Afc Mas. (14) 
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TABLE III: Relativistic deflection of light by some planetary satellites. 



Object 


Mass, 
iU g 


Radius, 
/v B , J^m 


Angular size, 
TZ- B , arcsec 


Grazing 
Pgr- , ^as 


1 fias critical radius 


dcrit, km 


^crit ) T^planet 


lo 


7.23 


1738 


0.570056 


25.48 


44 291 


0.63 ■7^J 


Europa 


4.7 


1620 


0.531353 


17.77 


28 793 


0.41 ■7^J 


Ganymede 


15.5 


1415 


0.464114 


67.11 


94 954 


1.34 ■7^J 


Callisto 


9.66 


2,450 


0.803589 


24.15 


59178 


0.84 ■ 7^J 


Rhea 


0.227 


675 


0.108468 


2.06 


1391 


0.02 • TLs 


Titan 


14.1 


2 475 


0.397715 


34.90 


86 378 


1.44 • lis 


Triton 


13 


1750 


0.082638 


45.51 


79 639 


3.17•7^Jv 



The necessary expression for dj^it is obtained from Eq. pOj) . Assuming that the angle xb is small and noting 
that rssinxs = d, we can write this equation as follows 9^^ ~ 2fiB/d, where fiB = 2GMb/c^ being relativistic 
gravitational radius of the body. As the effect of gravitational deflection light is inversely proportional to the impact 
parameter, then beginning from a certain value of the parameter, d^iti the deflection angle will be larger A^o; this 
value is given by the following expression: 

Different forms of the critical impact parameters d^^^^ for A6'o = 1 /ias are given in the Table HT] With the help of 
Eq. pS]). the results given in this table are easily scaled for any astrometric accuracy A6'o. 



2. Deflection of Light by Planetary Satellites 



One may expect that the planetary satellites will affect the astrometric studies if a light ray would pass in their 
vicinities. Just for completeness of our study we would like to present the estimates for the gravitational deflection of 
light by the planetary satellites and the small bodies in the solar system. The corresponding estimates for deflection 
angles, 6*^, and critical distances, dcrit are presented in the Table Hill Due, to the fact that the angular sizes for 
those bodies are much less than the smallest field of regard of the SIM instrument (e.g. FoR=l°), the results for the 
differential observations will be effectively insensitive to the size of the the two available FoRs. The obtained results 
demonstrate the fact that observations of these objects with that size of FoR will evidently have the effect from the 
relativistic bending of light. Thus, in Table Hill we have presented there only the angle for the absolute gravitational 
deflection in terms of quantities 0^. As a result, the major satellites of Jupiter, Saturn and Neptune should also be 
included in the model if the light ray passes close to these bodies. 



3. Gravitational Influence of Small Bodies 

Additionally, for the astrometric accuracy at the level of few yitas (i.g., A^o = Afc ^as), one needs to account for the 
post-Newtonian deflection of light due to rather a large number of small bodies in the solar system having a mean 
radius 




7^s > 624 J km. (16) 



The deflection angle for the largest asteroids Ceres, Pallas and Vesta for A/c = 1 are given in the Table HVl The 
quoted properties of the asteroids were taken from [2^, [s^l- Positions of these asteroids are known and they are 
incorporated in the JPL ephemerides. Other small bodies (e.g. asteroids, Kuiper belt objects, etc.) may produce a 
stochastic noise in the future astrometric observations with SIM; therefore, they should also be properly modeled. 



III. REGIONS WITH THE MOST GRAVITATIONALLY INTENSE ENVIRONMENTS FOR SIM 



The properties of the solar system's gravity field presented in the Tables U and |TT] suggest that the most intense 
gravitational environments in the solar system are those offered by the Sun and two planets, namely the Earth and 
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TABLE IV: Relativistic deflection of light by the asteroids. 



Object 


PB, g/cm^ 


Radius, km 


6^, fias 


Ceres 


2.3 


470 


1.3 


Pallas 


3.4 


269 


0.6 


Vesta 


3.6 


263 


0.6 


Class S 


2.1 ± 0.2 


TBD 


< 0.3 


Class C 


1.7 ± 0.5 


TED 


< 0.3 



TABLE V: Magnitudes of the gravitational deflection angle vs. the Sun-source separation angle xi©- 



Solar 
deflection 


smaU xi0, deg 


0".27 


0".5 


r 


2" 


5" 


10" 


15" 


6'^ , mas 


1728 


933.295 


466.639 


233.302 


93.271 


46.547 


30.932 


(56l« [15°], mas 


1698 


903.372 


437.663 


206.053 


70.176 


28.178 


15.734 


5eg [1°], mas 


1361 


622.212 


233.337 


77.787 


15.567 


4.254 


1.956 



Solar 
deflection 


large xiQ, deg 


20° 


40" 


45" 


50" 


60" 


70" 


80" 


90° 


6^ , mas 


23.095 


11.189 


9.832 


8.733 


7.053 


5.816 


4.853 


4.072 


56 [15°], mas 


10.180 


3.366 


2.778 


2.341 


1.746 


1.372 


1.122 


0.948 


56l« [1°], mas 


1.123 


0.297 


0.238 


0.195 


0.140 


0.107 


0.085 


0.071 



Jupiter. In this Section we will analyze these regions in more details. 



A. Gravitational Deflection of Light by the Sun 

From the expressions Eq. ([9]) and Eq. (jTS]) we obtain the relations for relativistic deflection of light by the solar 
gravitational monopole. The expression for the absolute astrometry takes the form: 

e = (7+l)^^ii^^= 4.072.i±i2^^ mas, (17) 

where xi© is the Sun-source separation angle, tau = 1 AU, and 7 = 1. Similarly, for differential astrometric 
observations one obtains: 

xaQ ^ Me sin ^(x20 - X10) , „„„ sin ^(x20 - X10) 

d9^^ = ("/ + l)— : — ; — = 4.072-^ — f — = mas, (18) 

TAu sm 2X10 ■ sm 5X20 sm 5X10 ' sm 5X20 

with Xi0iX20 being the Sun-source separation angles for the primary and the secondary stars correspondingly. We 
use two stars separated by the SIM's field of regard, namely X20 = X10 + 15 • The solar angular dimensions from the 
Earth' orbit are calculated to be TZq — 0°. 26656. This angle corresponds to a deflection of light to 1.75065 arcsec on 
the limb of the Sun. Results for the most interesting range of xi0 are given in the Table fVl 

Figure [4] shows a qualitative presentation of the solar gravitational deflection. The upper thick line on both plots 
represents the absolute astrometric measurements, while the other two are for the differential astrometry. Thus, the 
middle dashed line is for the observations over the maximal field of regard of the instrument FoR = 15°, the lower 
thick line is for FoR = 1°. 

One can also account for the post-post-Newtonian (post-PN) terms (e.g. cx G^) as well as the contributions due 
to other PPN parameters [l3|- Thus, in the weak gravity field approximation the total deflection angle 6'gr has an 
additional contribution due the post-post-Newtonian terms in the metric tensor. For the crude estimation purposes 
this effect could be given by the following expression [sj] : 

Ho.-™ - ^(7 + 1)^ j - 1 j [—^ ) . (19) 
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10 20 30 40 40 60 80 100 120 140 

Sun-source separation angle, xi0 [deg] Sun-source separation angle, xi© [deg] 



FIG. 4: Solar gravitational deflection of light. On all plots: the upper thick line is for the absolute astrometric measurements, 
while the other two are for the differential astrometry. Thus, the dashed line is for the observations over fleld of regard of 
FoR = 15°, the lower thick hue is for FoR = 1°. 



However, a quick look on the magnitudes of these terms for the solar system's bodies suggested that SIM astrometric 
data will be insensitive to the post-PN effects. The post-PN effects due to the Sun are the largest among those in the 
solar system. However, even for the absolute astrometry with the Sun-grazing rays the post-PN terms were estimated 
to be of order (5^^st-PN = 11 /^as. Note that the SIM solar avoidance angle is constraining the Sun-source separation 
angle as xi© > 45°. The post-PN effect is inversely proportional to the square of the impact parameter, thus reducing 
the effect to ^^^st-PN — 4.9 nanoarcseconds at the edge of the solar avoidance angle. Therefore, the post-PN effects 
will not be accessible with SIM. 



B. Gravitational Deflection of Light by Jupiter 

Astrometric measurement with SIM would have to account for the light bending by Jupiter [s^, One may 
obtain the expression, similar to Eq. p7p for the relativistic deflection of light by the Jovian gravitational monopole 
in the following form: 

e = (7 + l)^^ii^^ = 0.924944. i±^^ /.as, (20) 

with xij being Jupiter-source separation angle as seen by the interferometer at the distance rj from Jupiter. For the 
differential observations one will have expression, similar to that Eq. for the Sun: 

.nj / Mj sin^(x2,7 -Xij) nao/iQ/i/i ^^'^l(X2j - Xu) 

oOgr = (7 + 1)^ — — = 0.924944 • ^—^ — ^as, (21) 

O sm 2X1J • sm 2X2J sm2XiJ-sm2X2j 

where again xu? X2J are Jupiter-source separation angles for the primary and secondary stars correspondingly, X2J = 
Xij + j2 (and X2J = Xu + yfp for the narrow angle astrometry). The largest effect will come when SIM and Jupiter 
are at the closest distance from each other ^ 4.2 AU. Jupiter's angular dimensions from the Earth' orbit for this 
situation are calculated to be TZj = 23.24 arcsec, which correspond to a deflection angle of 16.419 mas. Results for 
some Xu are given in the Table IVll Note that for the light rays coming perpendicular to the ecliptic plane the Jovian 
deflection will be in the range: 5aij ^ (0.7 — 1.0) /xas! 

A qualitative behavior of the effect of the gravitational deflection of light by the Jovian gravity field is plotted in the 
Figure [51 As in the case of the solar deflection, the upper thick line on both plots represents the absolute astrometric 
measurements, while the other two are for the differential astrometry (the dashed line is for the observations over 
FoR = 15° and the lower thick hne is for FoR = 1°). 
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TABLE VI: Jovian gravitational monopole deflection vs. the Jupiter-source sky separation angle xu- 



Jovian 
deflection 


Jupiter-source separation angles xu, arcsec 


23.24" 


26" 


30" 


60" 


120" 


180" 


360" 


90" 


Sgr, mas 


16.419 


14.676 


12.719 


6.360 


3.180 


2.120 


1.060 


0.9 /las 


Se^^[15°], mas 


16.412 


14.669 


12.712 


6.352 


3.173 


2.113 


1.053 


0.2 ^as 


(56lg4[l°], mas 


16.313 


14.570 


12.614 


6.255 


3.077 


2.019 


0.964 


0.0 ^as 




C. Gravitational Deflection of Light by the Earth 

The deflection of light rays by the Earth's gravity field may also be of interest. The expressions, describing the 
relativistic deflection of light by the Earth' gravitational monopole are given below: 

e ^ + 1) g 1 + cosxie ^ . 1 + cosxie ^^^^ (22) 
^ & r® smxie smxi© 

with xi® being the Earth-source separation angle as seen by the interferometer at the distance from the Earth. 
Relation for the differential astrometric measurements was obtained in the form: 



s^a® ( , 1^G Me Sin2(x2©-Xie) rn^/is sm2(X2e-Xi©) .„„^ 

5C = (7 + 1)^ ^ — = 0.2446 • — //as, (23) 

re sm 2X1© • sm 3X2© sm 2X1© ' sm 3X2© 

where, as before, Xi©) X2© are the Earth-source separation angles for the primary and secondary stars correspondingly, 
X2© = Xi© + j2 (and X2© = Xiffi + ifo f^*^^ the narrow angle astrometry). The largest effect will come when SIM 
and the Earth are at the closest distance, say at the end of the first half of the first year of the mission, r® = 0.05 
AU. The Earth's angular dimensions being measured from the spacecraft from that distance are calculated to be 
"R?™ = 175.88401 arcsec, which correspond to a deflection angle of 573.75 /las. The summary of the deflection angles 
for sevral xi® are given in the Table IVTIl 

Figure [6] shows the expected variation in the magnitude of the Earth' gravity influence as mission progresses. The 
left plot presents results for the end of the first half of the year of the mission, when the spacecraft is at the distance 
of 0.05 AU from the Earth (the drift rate is 0.1 AU per year). The plot on the right side is for the end of the 5-th 
year of the mission, when SIM is at 0.5 AU from Earth. 
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TABLE VII: Solar relativistic deflection angle as a function of the Earth-source separation angle. 



SIM 
mission 


Xie, arcsec 


175.88 


200 


360 


r 


5" 


10" 


15" 


010, Atas 


573.8 


504.7 


280.3 


28.0 


5.6 


2.8 


1.9 


(5^1® [15°], pias 


571.9 


502.7 


278.5 


26.3 


4.2 


1.7 


1.0 


5^1® [1°], Atas 


547.0 


478.0 


254.8 


14.0 


0.9 


0.3 


0.1 



At the distance of 0.05 AU from the Earth At the distance of 0.5 AU from the Earth 




2 4 6 8 10 12 14 2 4 6 8 10 12 14 

Earth-source separation angle, xie [deg] Earth-source separation angle, xi® [deg] 

FIG. 6: Gravitational deflection of light in the proximity of the Earth. 



D. Constraints Derived EVom the Monopole Deflection of Light 

While analyzing the solar gravity field's influence on the future astrometric observations with SIM, we found 
several interesting situations, that may potentially put an additional navigational requirements. In this section we 
will consider these situations in a more detailed way. 

To measure gravitational deflection of light with an accuracy of A^o = Afc ^as, one needs to precisely determine 
the value of impact parameter of photon's trajectory with respect to the deflector. As before, we will present two 
types of necessary expressions, namely for absolute and differential observations. 

In the case of absolute astrometry we use expression for the deflection angle from Eq. (fTO|) and present it as 
^gr — (M-b/'^)(1 + cosxs), where again d = rg s'mxB and quantity iib = ^GMg/c^ being the relativistic gravitational 
radius of the body at question. One may ask a question - what uncertainty in the knowledge of the impact parameter, 
AcZ, will result in the astrometric error of A^o? The answer is given the following expression 

Ads ^A9,±. ^11^. (24) 
/is 1 + cos XlB 

The corresponding result for differential observations may be obtained with the help of Eq. as: 



1 + tan ^ • cot i(x2B - Xib) ■ (25) 



2^lB 

Similarly, the uncertainty in determining the barycentric distance is obtained from Eq. (jlOp leading to expression: 

A.. = A^o^- . (26) 

/is 1 + COSXlB 

Note that, when differential observations are concerned, uncertainty in barycentric position Ar^^*' does not produce 
new constraints significantly different from those derived from Eq. (|26p . Looking at the results presented in the 
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TABLE VIII: Required accuracy of barycentric positions and impact parameters for astrometric observations with accuracy of 
1 /las. The Earth is taken at the distance of 0.05 AU from the spacecraft; moon's position accuracy is for the geocentric frame. 



Solar 


Required knowledg 


e: grazing rays 


Required knowledge: 


differential astrometry 


system's 


Distance, 


Impact parameter 


Impact param. [15"] 


Impact param. [1"] 


object 


(Jr a , km 


(Jdg , km 


, mas 


(Jrg , km 


, mas 


cTrg , km 


, mas 


Sun 


85.45 


0.39 


0.55 


0.40 


0.55 


0.50 


0.69 


Sun at 45° 


1.5x10* 


7.6x 10^ 


10". 49 


3.81 X 10 


52". 53 


4.45x 10 


613 .66 


Moon 


2.8x 10 


({I 1 A 


1 " sc. 

i .00 


f<7 1 7 
i .1 ( 


1 .OO 


Do.UU 


1 " 

1 .OO 


Mercury 


i.iX iU 


29.39 


66.16 


29.41 


66.16 


29.45 


66.25 


Venus 


8.4x10* 


12.18 


61.00 


12.28 


61.20 


12.38 


62.70 


Earth-Moon 


1.3x10* 


11.14 


306.55 


11.15 


307.47 


11.66 


321.54 


Mars 


6.8x10^ 


29.29 


77.11 


29.30 


77.14 


29.37 


77.33 


Jupiter 


3.8x10* 


4.31 


1.42 


4.32 


1.42 


4.34 


1.42 


Jupiter at 30" 


4.9x10* 


7.14 


2.34 


7.20 


2.36 


7.25 


2.38 


Saturn 


2.2x10^ 


10.32 


1.66 


10.34 


1.66 


10.36 


1.66 


Uranus 


1.2x10® 


11.27 


0.86 


11.28 


0.86 


11.29 


0.86 


Neptune 


1.7x10® 


10.04 


0.47 


10.04 


0.47 


10.04 


0.47 


Pluto 


2.0x10^* 


1 133.92 


40.7 


1 133.93 


40.7 


1133.96 


40.7 



TABLE IX: The best known accuracies of barycentric positions and masses for the solar system's objects [3^. l35l, [s^ . 



Solar 
system's 
object 


Knowledge of barycentric position 


Knowledge of 
planetary masses, 
AMb/Mb 


Best known, 


Method used 
for determination 


, km 


(Trg , mas 


Sun 


362/725 


0".5/l".0 


Optical meridian transits 


3.77x10"'" 


Moon 


27 cm 


7.4 fias 


LLR, 1995 


1.02x10"® 


Mercury 


1 


2.25 


Radar ranging 


4.13x10"^^ 


Venus 


1 


4.98 


Radar ranging 


1.23x10"'^ 


Earth 


1 


27.58 


Radar ranging 


TBDxlO"® 


Mars 


1 


2.63 


Radar ranging 


2.33x10"® 


Jupiter 


30 


9.84 


Radar ranging 


7.89x10"^^ 


Saturn 


350 


56.24 


Optical astrometry 


2.64x10"® 


Uranus 


750 


57.00 


Optical astrometry 


3.97x10"® 


Neptune 


3 000 


141.67 


Optical astrometry 


2.19x10"® 


Pluto/Charon 


20 000 


717.40 


Photographic astrometry 


0.014 



Table I Villi one may see that for an accuracy of A6'o — 1 /ias our estimates require the knowledge of the solar impact 
parameter with the accuracy of ~ 0.4 km (grazing rays), that for Jupiter with the accuracy of ~ 4 km and other 
big planets with the accuracy of about 10 km. Table IIXI shows a comparison of these derived requirements on the 
barycentric positions of the solar system's bodies with the accuracy of their current determination. 

One may see that the present accuracy of knowledge of the inner planets' positions from the Table HXl is given by 
the radio observations and it is even better than the level of relativity requirements given in the Table IVTIIl However, 
the positional accuracy for the outer planets is below the required level. The SIM observation program should include 
the astrometric studies of the outer planets in order to minimize the errors in their positional accuracy determination. 
Thus, in order to get the radial uncertainty in Pluto's ephemeris with accuracy below 1000 km, it is necessary only 4 
measurements of Pluto's position, taken sometime within a week of the stationary points, spread over 3 years. Each 
measurement could be taken with an accuracy of about 200 /las, as suggested by [331 • Additionally, one will have to 
significantly lean on the radio observations in order to conduct the reduction of the optical data with an accuracy of 
a few ^as. For this reason one will have to use t he p recise catalog of the radio-sources and to study the problem of 
the radio and optical reference frame ties [l^, 113, [H, HI] . 

The estimates, presented here were given for static gravitational field. Analysis of a real experimental situation 
should consider a non-static gravitational environment of the solar system and should include the description of light 
propagation in a different reference frames involved in the experiment [2^ . Additionally, the observations will be 
affected by the relativistic orbital dynamics of the spacecraft [2l| . 
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TABLE X: Higher gravitational coefficients for solar system bodies (|http: / /nssdc .gsf c .nasa.gov/plaiietary/f actsheet/[ |. 



Solar system's 


■rH 


J4 , 




object 


X 10~^ 


X 10"*^ 


X 10"" 


Sun 


0.17 ±0.017 


— 


— 


Sun at 45° 


— 


— 


— 


Moon 


202.2 


-0.1 


— 


Mercury 


60. 


— 


— 


Venus 


4.5 


-2.1 




Earth 


1082.6 


-1.6 


0.5 


Mars 


1 960.45 






Jupiter 


14 738±1 


-587±5 


34±50 


Saturn 


16 298±50 


-915±80 


103.0 


Uranus 


3 343.43 






Neptune 


3411. 






Pluto 









IV. DEFLECTION OF LIGHT BY HIGHER GRAVITY MULTIPOLES 



In order to carry out a complete analysis of the relativistic light deflection one should account for other possible 
terms in the expansion ([T]) that may potentially contribute to this effect. These terms are due to non-sphericity 
and non-staticity of the body's gravity field [5|]. Here we will consider several of them, namely those due to higher 
gravitational multipoles of the celestial bodies. 



A. Gravitational quadrupole deflection of light 

Although a complete three-dimensional deflection of light must be considered for a real experiment, for the purposes 
of this paper, we consider only two dimensional case. In this case, quadrupole term may be given as P. [2ll. |40| : 

I, , ,,4GMs , 2 ^ 2\ fT^B 



where J2 is the second zonal harmonic of the body under question, s = [sx, Sy,Sz) is the unit vector in the direction 
of the light ray propagation and vector d — d{ux, Uy,Uz) is the impact parameter. A similar expression may also be 
obtained for differential observations. For estimation purposes, this formula may be given as follows (d — tb suyxb)'- 



„ 4:HbJ2bTZ% 



sin^ XiB sin^ X2B 



(28) 



The corresponding effects for the deflection of light by the quadrupole mass moments within the planets of the solar 
system are given in the Table I XII The effect depends on a number of different instantaneous geometric parameters 
defining the mutual orientation of the vector of the light propagation, position of the planet in orbit, the orientation 
of the axes defining it's figure, etc. A mission-independent modeling in the static gravitational regime has been done 
[s^l; the modeling of the dynamic regime was presented in An effort to develop a SIM-specfic model is well 
justified. 

The quadruple deflection of light depends on the third power of the inverse impact parameter with respect to the 
deflecting body, Eq. ([77|) . SIM will measure this effect directly for many celestial bodies. At the expected level of 
accuracy the knowledge of jovian atmosphere, the magnetic field fluctuations, etc., may contribute to the errors in 
the experiment [23| . A detailed study of these effects is given in [4l[ . 

Table |Xl] presents estimates of the magnitudes of the relativistic deflection of light by the planetary quadrupole mass 
moments for solar system bodies. Based on the sizes of these effects, one would have to account for the quadrupole 
component of the gravity fields when observations will be conducted in the vicinity of the outer planets. In addition, 
the infiuence of the higher harmonic may be also of interest. Thus, Table IXII shows the estimates of some higher 
gravitational multipole moments of Jupiter and Saturn. We will discuss the deflection by the J2 and J4 coefficients 
of the jovian gravity field in terms of Jupiter-source separation angle xu- An expression, similar to that of Eq. (^0]) 
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TABLE XI: Relativistic deflection of light by the planetary quadrupole mass moments for solar system bodies. 



OOldi byoLclli 

object 


X 10-** 


f J2> 

/ias 




A/9« [1 c:°l 
"t^Ja 1-'-° J' 

/ias 


"^.121^ J' 

/ias 


Sun 


0.17 ± 0.017 


0.3 




0.3 


0.3 


Moon 


202.2 


2x10"^ 


— 


2x10-2 


2 x 10-^ 


Mercury 


60. 


5x10"^ 








Venus 


4.5 


2x10"^ 


— 


— 


— 


Earth 


1082.6 


0.6 




0.6 


0.6 


Mars 


1 960.45 


0.2 


— 


0.2 


0.2 


Jupiter 


14 738±1 


242.0 


98". 12 - 144". 81 
6.23 TZj 


242.0 


242.0 


Saturn 


16 298± 50 


94.6 


35" .62 - 43". 93 
4.56 7^s 


94.6 


94.6 


Uranus 


3 343.43 


7.3 


3". 25 - 3". 61 
1.94 Tlu 


7.3 


7.3 


Neptune 


3411. 


8.5 


2". 23 - 2".42 
2.04 Un 


8.5 


8.5 


Pluto 













for the monopole deflection, may be given as: 

fff^ = 3.46058 x 10-1° — I ^las. (29) 

sm'" Xij 

Jupiter's angular dimensions from the Earth are calculated to be TZj = 23.24 arcsec, which correspond to a deflection 
angle of 242 /ias. The deflection on the multipolcs for some Xij is given in the Table IXlIl 



TABLE XIL Deflection of light by the Jovian higher gravitational coefficients. 



Jovian 
deflection 


Xu, arcsec 


23". 24 


26" 


30" 


35" 


40" 


50" 


120" 


61;},, /ias 


242 


173 


112 


71 


47 


24 


1.8 


<56»,jJ15°], /ias 


242 


173 


112 


71 


47 


24 


1.8 


Oj^, /ias 


9.6 


5.5 


2.7 


1.3 


0.6 


0.2 


0.0 



TABLE XIIL Deflection of light by the Saturnian higher gravitational coefficients. 



Saturnian 
deflection 


Xis, arcsec 


9". 64 


12" 


15" 


20" 


25" 


30" 


35" 


e%, /ias 


94.7 


49.1 


25.1 


10.6 


5.4 


3.1 


2 


5e5J15°], /ias 


94.7 


49.1 


25.1 


10.6 


5.4 


3.1 


2 


Mas 


5.3 


1.8 


0.6 


0.1 









Similar studies are important for Saturn (see [J] for details). In terms of the Saturn-source separation angle xis 
the saturnian quadrupole deflection mat be estimated with the help of the following expression: 

OT" = 9.66338 X 10"^^ — ^ ^as. (30) 

sm-" xis 

The Saturn's angular dimensions from the Earth' orbit are calculated to be TZs = 9.64 arcsec, which correspond to a 
deflection angle of 94.7 /ias. The corresponding estimates for the deflection angles are presented in the Table IXlIIl 

As a result, for astronomical observations with accuracy of about 1 /ias, one will have to account for the quadrupole 
gravitational fields of the Sun, Jupiter, Saturn, Neptune, and Uranus. In addition, the influence of the higher 
harmonics may be of interest. For example some of the moments for Jupiter and Saturn are given in the Table 1x1 
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Higher multipoles may also influence the astrometric observations taken close to these planets. Thus, for both 
Jupiter and Saturn the rays, grazing their surface, will be deflected by the fourth zonal harmonic J4 as follows: 
59 « 9.6 /ias, 56^^ ~ 5.3 /ias. In addition, the contribution of the Jq for Jupiter and Saturn will deflect the grazing 
rays on the angles ~ 0.8 /las, 59^j^ ~ 0.6 /J,as. The contribution of J4 is decreasing with the distance from the 
body as and contribution of Je as As a result the deflection angle will be less then 1 /las when d > 1.6 TZ, 
where TZ is the radius of the planet. 

Using Eq. ([?f|) . one can derive the expression for the critical distance rf}^^^ for the astrometric observations in the 
regime of quadrupole deflection of light with accuracy of A^o — /las. Indeed, approximating this equation as 
~ {AfiB/d^)J2Ti'^ and solving it for d one obtaines the following result: 



d7r = nB 



-Rb A9o 



(31) 



The critical distances for the relativistic quadrupole deflection of light by the solar system's bodies for the case of 
Afc = 1 presented in the Table [ 



B. Gravito-Magnetic Deflection of Light 



Besides the gravitational deflection of light by the monopole and the quadrupole components of the static gravity 
fleld of the bodies, the light ray trajectories will also be affected by the non-static contributions from this field. It is 
easy to demonstrate that a rotational motion of a gravitating body contributes to the total curvature of the space-time 
generated by this same body. This contribution produces an additional deflection of light rays 0, |43| on the 
angle 

S9^-=^{j + l)^S{s-d), (32) 

where S is the body's angular momentum. 

The most signiflcant contributions of gravito-magnetic deflection of light by the bodies of the solar system arc the 
following ones: i) the solar defiection amounts to 69'^ = ±(0.7 — 1.3)/ias (the first term listed is for a uniformly 

rotating Sun; the second is for the Dicke's model Hi); ii) jovian rotation contributes 69'^ = ±0.2 /xas; and iii) 
Saturn's rotation S9'i°' = ±0.04 fias. Thus, depending on the model for the solar interior, solar rotation may produce 

o 

a noticeable contribution for the grazing rays. The estimates of magnitude of deflection of light ray's trajectory, 
caused by the rotation of gravitating bodies demonstrate that for precision of observations of 1 /ias it is sufficient to 
account for influence of the Sun and Jupiter only. 

The relativistic gravito-magnetic deflection of light has never been directly tested before. However, because of the 
smallness of the magnitudes of corresponding effects in the solar system and SIM's operational mode that limits the 
viewing angle for a sources as X10 > 45°, SIM will not be sensitive to this effect. 



V. ASTROPHYSICS INVESTIGATIONS WITH SIM 



A. Astrometric Test of General Relativity 



The Eddington parameter 7 in Eq. ^ , whose value in general relativity is unity, is perhaps the most fundamental 
PPN parameter [2J], in that ^(1—7) is a measure, for example, of the fractional strength of the scalar gravity 
interaction in scalar-tensor theories of gravity [1^. Currently, the most precise value for this parameter, 7 — 1 = 
(2.1 ± 2.3) X 10~^, was obtained using radio-metric tracking data received from the Cassini spacecraft [25|| during 
a solar conjunction experiment. This accuracy approaches the region where multiple tensor-scalar gravity models, 
consistent with the recent cosmological observations [i^ , predict a lower bound for the present value of this parameter 
at the level of (1 — 7) ^ 10~^ — 10~^ j45j . Therefore, improving the measurement of this parameter would provide the 
crucial information separating modern scalar-tensor theories of gravity fromgeneral relativity, probe possible ways 
for gravity quantization, and test modern theories of cosmological evolution [l|,[3]. 

The reasons above led to a number of speciflc space experiments dedicated to measurement of the parameter 7 
with a precision better than 10~^ to 10~^ [H, 0- Note that SIM will operate at this level of accuracy and, therefore, 
the Eddington's parameter 7 will have to be included into the future SIM's astrometric model and the corresponding 
data analysis. 



17 

1. Solar Gravity Field as a Deflector 

To model the astrometric data to the nominal measurement accuracy will require including the effect of general 
relativity on the propagation of light. In the PPN framework, the parameter 7 would be part of this model and could 
be estimated in global solutions. The astrometric residuals may be tested for any discrepancies with the prescriptions 
of general relativity. To address this problem in a more detailed way, one will have to use the astrometric model 
for the instrument including the information about it's position in the solar system, it's attitude orientation in the 
proper reference frame, the time history of different pointings and their durations, etc. This information then should 
be folded into the parameter estimation program that will use a model based on the expression, similar to that given 
by Eq. (ini). 

To estimate the expected accuracy of the parameter 7, we use Eq. ()13p to assume that the single astrometric 
measurement may be able to determine this parameter with accuracy: 

A7 = A^„^^^^i^^^^^, (33) 
Ai0 sm2(X20-Xi0) 

where A^o is the astrometric error of the measurement. 

The relativity test will be enhanced by scheduling measurements of stars as close to the Sun as possible. Although 
SIM will never be able observe closer to the Sun than 45°, it will allow for an accurate determination of this PPN 
parameter. For the accuracy of A9a = 1 fias at the rim of the solar avoidance angle of xiq = 45°, one could determine 
this parameter with an accuracy ~ 7.19 x 10~* in a single measurement. Assuming Gaussian error distribution, the 
accuracy of this experiment will improve as 1/Vn, where N is the number of independent observations. Therefore, by 
performing differential astrometric measurements with an accuracy of AOq = 1 /ias over the instrument's FoR=15°, 
at the end of the mission (with N ^ 10000 observations) SIM may reach the accuracy of cr^ ~ 7.2 x 10^^ in 
astrometric test of general relativity in the solar gravity field. SIM will provide this precision as a by-product of its 
astrometric program, thus allowing for a factor of 3 improvement of the currently best Cassini's 2003 result [1^. Such 
a measurement improves the accuracy of the search for cosmologically relevant scalar-tensor theories of gravity by 
looking for a remnant scalar field in today's solar system. 



2. GR Test in the Jovian and Earth' Gravity Fields 



One could also perform a relativity experiment with Jupiter and the Earth. In fact, for the proposed SIM's observing 
mode, the accuracy of determining of the parameter 7 may be even better than that achievable with the Sun. Indeed, 
with the same assumptions as above, one may achieve a single measurement the accuracy of I7 — 1| ^ 4.0 x 10""* 
determined via deflection of light by Jupiter. As the astrometric observations in the Jupiter's vicinity will require 
careful planning thereby minimizing the number of possible independent observations. As a result, the PPN parameter 
7 may be obtained with accuracy of about cr^ ~ 1.3 x 10~^ with astrometric experiments in Jupiter's gravity field 
(note that only N ^ 1000 needed). 

Lastly, let us mention that the experiments conducted in the Earth's gravity field, could also determine this 
parameter to an accuracy I7 — 1| ~ 8.9 x 10""^ in a single measurement (which in return extends the measurement 
of the gravitational bending of light to a different mass and distances scale, as shown by [l3|)- One may expect a 
large statistics gained from both the astrometric observations and the telecommunications with the spacecraft. This, 
in return, will significantly enhance the overall solution for 7 obtained in the Earth' gravitational environment. 



B. Solar Acceleration Towards the Galactic Center 



The Sun's absolute velocity with respect to a cosmological reference frame was measured photometrically: it shown 
up as the dipole anisotropy of the cosmic microwave background j46| . The Sun's absolute acceleration with respect 
to galactic frame can be measured astrometrically: it will show up as a dipole vector harmonic in the global pattern 
of proper motion of quasars. 

The aberration due to the solar system's galactocentric motion will not be observable because its main contribution 
is static. However, the rate of this aberration will produce an apparent proper motion for the observed sources p^ . 
Indeed, the solar system's orbital velocity around the galactic center causes an aberrational affect of the order of 2.5 
arcmin. All measured star and quasar positions are shifted towards the point on the sky having galactic coordinates 
I = 90°, b = 0°. For an arbitrary point on the sky the size of the effect is 2.5 sin 77 arcmin, where rj is the angular 
distance to the point / = 90°, 6 = 0°. The acceleration of the solar system towards the galactic center causes this 
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aberrational effect to change slowly. This leads to a slow change of the apparent position of distant celestial objects, 
i.e. to an apparent proper motion. 

Let us assume a solar velocity of 220 km/sec and a distance of 8.5 kpc to the galactic center. The orbital period 
of the Sun is then 250 million years, and the galactocentric acceleration takes a value of about 1.75 x 10"^'^ km/sec^. 
Expressed in a more useful units it is 5.5 mm/s/yr. A change in velocity by 5.5 mm/sec causes a change in aberration 
of the order of 4 /las. The apparent proper motion of a celestial object caused by this effect always points towards 
the direction of the galactic center. Its size is 4 sin 77 /ias/yr, where i] is now the angular distance between the object 
and the galactic center. 

The above holds in principle for quasars, for which it can be assumed that the intrinsic proper motions (i.e. those 
caused by real transverse motions) are negligible. A proper motion of 4 /las/yr corresponds to a transverse velocity 
of 2 X 10^ km/sec at z 0.3 for i?o=100 km/sec/Mpc, and to 4 x 10^ km/sec for Hq = 50 km/sec/Mpc. Thus, all 
quasars will exhibit a distance- independent steering motion towards the galactic center. Within the Galaxy, on the 
other hand, the effect is drowned in the local kinematics: at 10 pc it corresponds to only 200 m/sec. 

However, for a differential astrometry with SIM this effect will have to be scaled down to account for the size of 
the field of regard ^2)], namely 2sin^ — 2sin^ = 0.261. This fact is reducing the total effect of the galactocentric 
acceleration to only 1 sin rj /las/yr and, thus, it makes the detection of the solar system's galactocentric acceleration 
with SIM to be a quite problematic issue. 

Discussion 

General relativistic deflection of light produces a significant contribution to the future astrometric observations 
with accuracy of about a few /ias. In this paper we addressed the problem of light propagation on the gravitational 
field of the solar system. It was shown that for high accuracy observations it is necessary to correct for the post- 
Newtonian deflection of light by the monopole components of gravitational fields of a large number of celestial bodies 
in the solar system, namely the Sun and the nine planets, together with the planetary satellites and the largest 
asteroids (important only if observations are conducted in their close proximity) . The most important fact is that the 
gravitational presence of the Sun, Jupiter and the Earth should be always taken into account, independently on the 
positions of these bodies relative to the interferometer. It is worth noting that the post-post Newtonian effects due 
to the solar gravity are unlikely to be accessible with SIM. This effect as well as the effect of gravitational deflection 
of light caused by the mass quadrupole term of the Sun are negligible at the level of expected accuracy. However, 
deflection of light by some planetary quadrupoles may have a big impact on the astrometric accuracy. Thus, the 
higher gravitational multipoles should be taken into account when observations are conducted in the close proximity 
of two bodies of the solar system, notably Jupiter and Saturn. 

We emphasized the need of development of a general relativistic model for SIM observables to enable the mission 
to improve the current astrometric accuracy by a factor of over 1000. This model would have to account for a 
number of dynamical effects both external to the spacecraft (e.g., motion with respect to the solar system barycentric 
reference frame, effects of time- varying gravitational field in the solar system (due to planetary motion nd rotation) 
on light propagation, various interplanetary media effects, etc.) and internal to the spacecraft (e.g., systematic effects 
introduced by the spacecraft itself). Some of this work has already begun in the context of the development relativistic 
reference frames for the need of future high- precision observations However, a lot more efforts is needed; this 
paper intends to motivate initiation of such a work in the near future. 
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